Entanglement entropy of two disjoint intervals in 
conformal field theory 



Pasquale Calabrese^, John Cardy^, and Erik Tonni^ 

^ Dipartimento di Fisica deU'Universita di Pisa and INFN, Pisa, Italy. 

^ Oxford University, Rudolf Peierls Centre for Theoretical Physics, 1 Keble 

Road, Oxford, OXl 3NP, United Kingdom and All Souls College, Oxford. 

Abstract. 

We study the entanglement of two disjoint intervals in the conformal field 
theory of the Luttinger liquid (free compactified boson). Trp^ for any integer n 
is calculated as the four-point function of a particular type of twist fields and the 
final result is expressed in a compact form in terms of the Riemann-Siegel theta 
functions. In the decompactification limit we provide the analytic continuation 
valid for all model parameters and from this we extract the entanglement entropy. 
These predictions are checked against existing numerical data. 
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1. Introduction 

The interest in quantifying the entanglement in extended quantum systems has been 
growing in recent times at an impressive rate, mainly because of its ability in detecting 
the scaling behaviour in proximity of quantum critical points (see e.g. Refs. [Tl[51|31ll] 
as reviews) . A particularly useful measure of entanglement in the ground-state of an 
extended quantum system is the entanglement entropy Sa- It is defined as follows. Let 
p be the density matrix of a system, which we take to be in the pure quantum state |^) , 
p= |\|/)(\1>|. Let the Hilbert space be written as a direct product H = Ha ®'Hb- A's 
reduced density matrix is pA — Tr^ p- The entanglement entropy is the corresponding 
Von Neumann entropy 



and analogously for Sb- When p corresponds to a pure quantum state Sa = Sb- 

The entanglement entropy revealed to be an optimal indicator of the critical 
properties of an extended quantum system when A and B correspond to a spatial 
bipartition of the systems. Well-known and fundamental examples are critical one- 
dimensional systems in the case when A is an interval of length i in an infinite line. 
In this case, the entanglement entropy follows the scaling [SJ [SI [7] 



where c is the central charge of the associated conformal field theory (CFT) and c'^ 
a non-universal constant. Away from the critical point, Sa saturates to a constant 
value [5] proportional to the logarithm of the correlation length |7]. This scaling 
allows to locate the position (where Sa diverges by increasing I) and a main feature 
(by the value of the central charge c) of quantum critical points displaying conformal 
invariance. 

The central charge is an ubiquitous and fundamental characteristic of the CFT, 
but it docs not always identify unequivocally the universality class of the transition. 
A relevant class of relativistic quantum field theories are the so-called Luttinger 
liquids, which describe an enormous number of physical systems of experimental 
and theoretical interest. Just to quote a few, the one-dimensional Bose gases with 
repulsive interaction, the (anisotropic) Heisenberg spin chains, carbon nanotubes are 
all described by Luttinger liquid theory in their gapless phases. Via bosonization, all 
these models can be written as free bosonic field theories with c = 1. The different 
universality classes are distinguished by the compactification radius R of the bosonic 
field, that corresponds to experimentally measurable critical exponents. 

The entanglement entropy of a single block of length according to Eq. ^ is 
transparent to the value of the compactification radius, because it depends only on 
c. In Ref. [5] it has been shown that instead the entanglement entropy of disjoint 
intervals depends explicitely on i?, and so it encodes universal properties of the CFT 
that are hidden in the entanglement of a single block. (Oppositely in 2D systems 
with conformal invariant wave-function, the entanglement entropy of a single region 
depends on R [10].) Eq. ([2]) in a CFT is calculated by a modification of the replica 
trick of disordered systems [5lf7]. In fact, one first calculates Trp^ for integral n, that 
results to be 



Sa 



Tr PA log PA , 



(1) 



Sa = 7:^og- + d- 



(2) 



a 




(3) 
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that is easily analytically continued to any complex value of n, and then Sa = 
— lim„^i 9„Trp^ gives Eq. The reason of this way of proceeding is that for 

integral n, Trp^ is the partition function on an n-sheeted Riemann surface obtained 
by joining consecutively the n sheets along region A (see next section for details). 
We will refer to this surface as Tin,N, where N is the number of disjoint intervals 
composing A. {TZn,N are fully defined by the 2N branch points Uj and Vj). In the 
case of a single interval, TZn,i is easily uniformised to the complex plane by a simple 
conformal mapping. Then the powerful tools of CFT give Eq. ([3]). 

However, when the subsystem A consists of several disjoint intervals, the analysis 
becomes more complicated. In Ref. [7] (and also in [8]), based on a uniformising 
transformation mapping TZn,N into the complex plane, a general result for Trp^ has 
been given. However, this result is in general incorrect. In fact, the surface TZn,N has 
genus (n — 1)(-/V — 1) and so for ^ 1 cannot be uniformised to the complex plane 
(at the level of the transformation itself, this has been discussed in some details [TT]). 
The case n = N = 2 has the topology of a torus, whose partition function depends on 
the whole operator content of the theory and not only on the central charge (see e.g. 
[12]). Consequently the simple formulas of Ref. [7] cannot be generally correct. The 
partition functions on Riemann surfaces with higher genus are even more complicated. 

We consider here the case of two disjoint intervals A = AiU A2 = [ui, vi] U [u2, V2] 
defining the surface TZn,2- By global conformal invariance Trp^ can be written as 

Trpl . . el (-. :) ^„(x) , (4) 

\\U1 - Vi\\u2 - V2\\U1 - V2\\U2 - Vi\ J 

where x is the four-point ratio (for real Uj and Vj , x is real) 
^ _ jui - vi){u2 - V2) 

(Ui - U2)ivi - V2) ' 

This can be written as 

^7J„,2 = Z^„_2^n{x) , (6) 

where Z^^^ is the incorrect result in Ref. [7j. We normalised such that TniO) = 1. 
The function !Fn{x) depends explictly on the full operator content of the theory and 
must be calculated case by case. 

In Ref. [5], using old results of CFT on orbifolded space [131 E], T2{x) has been 
calculated for a free boson compactified on a circle of radius R 

e^{r,T)0,{T/rj) 

where r is pure-imaginary, and is related to x via x = [02(t)/03(t)]'*. 9i, are Jacobi 
thcta functions. 77 is a universal critical exponent proportional to the square of the 
compactification radius R (in Luttinger liquid literature rj = 1/{2K)). 
The main result of this paper is !Fn{x) for generic integral n > 1: 

_ 9(ohr)9(o|r/,) 
[e(o|r)]2 ' 

where F is an (n — 1) x (n — 1) matrix with elements 

^ (9) 



fc = i 



2z "^ . ( k\ 

— ^ ^ Sm I TT- 1 Pk/nCOS 



2TT—{r — s) 
n 
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and 

= Fy{x) = 2Fi{y,l-y;l;x). (10) 

Tj is exactly the same as above, while 9 is the Riemann-Siegel theta function 

e(z|r) ^ 

exp [iTrm' • r • TO + 27riTO* • z] , (11) 

with z a generic complex vector with n — 1 components. O(0|r) for n — 1 = 1 reduces 
to the Jacobi 9z{t = i(3i/2), and so Eq. ([8]) reproduces Eq. ([7]). In the above, 
and hereafter, the dot (•) denotes the matrix product and the superscript t (*) the 
transposition. Eq. ([5]) is manifestly invariant under rj l/rj, as numerically observed 
[9]. It is also invariant under x — *■ 1 — x (even if not manifest in this form). For 77 = 1 
we have !Fn{x) = 1 and the result ^ in [7] is then correct. This equality carries 
over to the analytic continuation and then to the entanglement entropy, confirming 
what observed numerically [9]. 

Unfortunately we have been not yet able to analytically continue this result to real 
n for general values of rj and x, and so to obtain the entanglement entropy. However 
we managed to give asymptotic expressions for small and large rj that compare well 
with numerics. 

It is worth to recall that in the case of two intervals, the entanglement entropy 
measures only the entanglement of the two intervals with the rest of the system. It 
is not a measure of the entanglement of one interval with respect to the other, that 
instead requires the introduction of more complicated quantities because Ai U A2 is 
in a mixed state (see e.g. Refs. [Elll^ for a discussion of this and examples). 

The paper is organized as follows. In Sec. [2] we recall how to obtain the 
entanglement entropy within CFT on Riemann surfaces and the usefulness of twist 
fields. In Sec. [3]we introduce the free compactified boson and fix all our notation. The 
following section [4] is the core of paper where Eq. (j4]) is derived. This section requires 
a good knowledge of CFT on orbifolds (to make the paper self-contained we also have 
a long I Appendix A] on this, where the results of Ref. [13] that we used are explained). 
The reader uninterested in the derivation can skip this section, to read directly Sec. 
[S] where a partial analytic continuation of Eq. ^ is performed and the consequences 
are discussed with particular attention to the comparison with the numerical results 
in Ref. [9j. Several appendices contain most of the technical parts of the paper. 



2. Entanglement entropy and Riemann surfaces 

Given a quantum field theory whose dynamics is described by the Hamiltonian H, the 
density matrix p in a thermal state at inverse temperature /3 may be written as a path 
integral in the imaginary time interval (0, (3) 

piiMlWA) = J my, r)] J] 5{cb{y, 0) -</.;,) H ^('^(j^' ^) " '^-) ' (^2) 

X X 

where Z{f3) — TTe~^^ is the partition function, the euclidean action is Se = Jq Ldr, 
with L the euclidean lagrangian. Here the rows and columns of the density matrix are 
labelled by the values of the fields at r = 0, (3. 

The normalisation factor of the partition function ensures that Trp = 1, and 
is found by setting {<j>x} — {4>'x} ^^'^ integrating over these variables. In the path 
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integral, this has the effect of sewing together the edges along t — and t — P to 
form a cylinder of circumference (3. Now let A be a subsystem consisting of the points 
X in the disjoint intervals {ui,vi), . . . , {un,vn)- An expression for the the reduced 
density matrix pA is obtained from (jl2p by sewing together only those points x which 
are not in A. This has the effect of leaving open cuts, one for each interval {uj,Vj), 
along the line t = 0. 

We may then compute Tr p\, for any positive integer n, by making n copies of the 
above, labelled by an integer j with 1 < j < n, and sewing them together cyclically 
along the the cuts so that 4>j(x,T — P") — <pj^i{x,T — 0^) and 4>n{x,T = j3~) = 
cj)i{x,T = 0+) for all x G A. This defines an n-sheeted Riemann surface depicted for 
n = 3 and in the case when A is formed by two disjoint intervals in Fig. [TJ The 
partition function on this surface will be denoted by Zn{A) and so 

TrPA = ^. (13) 

When the right hand side of the above equation has a unique analytic continuation to 
Ren > 1, its first derivative at n = 1 gives the required entropy 

Notice that even before taking the replica limit, these partition functions give the 
Renyi entropies 

5i"^ = Y^logTrpl. (15) 

Since the lagrangian density does not depend explicitly on the Riemann surface 
TZn.N as a consequence of its locality, it is expected that the partition function can 
be expressed as an object calculated from a model on the complex plane C, where 
the structure of the Riemann surface is implemented through appropriate boundary 
conditions around the points with non-zero curvature. Consider for instance the simple 
Riemann surface TZn,i needed for the calculation of the entanglement entropy of a 
single interval [uiiWi], made of n sheets sequentially joined to each other on the 
segment x S [ui,vi], t — 0. We expect that the associated partition function in a 
theory defined on the complex plane z — x -\- ir can be written in terms of certain 
"fields" at z = vi and z — ui. 

The partition function (here C[(p]{z, z) is the local lagrangian density) 



Ti„,N = / M'i5]7?,„,„ exp - / dzdz£[ip]{z,z) 



(16) 
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essentially defines these fields (i.e. it gives their correlation functions, up to a 
normalisation independent of their positions). In order to work with local fields (see 
for an extensive discussion [TT]), it is useful to move the intricate topology of the 
world-sheet (i.e. the space where the coordinates x,t lie) TZn,N to the target space 
(i.e. the space where the fields lie). Let us consider a model formed by n independent 
copies of the original model. The partition function (fT6|) can be re- written as the path 
integral on the complex plane 



Zn„ « = / [dfi ■ ■ ■ dLpn]c exp - / 



dzdz {C[(pi]{z, z) + . . . + C[ipn] [z, z)) 



(17) 



where with we indicated the restricted the path integral with conditions 

N 

Cuj,vj-- (y5i(a;,0+) = (^i+i(x,0") , x e [J w^], z = l,...,n, (18) 
where we identify n + i = i. The lagrangian density of the multi-copy model is 

£^"^[V3l, . . . , ipn][x,T) = C[ipi]{x,T) + ... + C[ipn]{x,T), 

SO that the energy density is the sum of the energy densities of the n individual copies. 
Hence the expression p7|) does indeed define local fields at (ui,0) and (wi,0) in the 
multi-copy model [T7] . 

The local fields defined in (|17p are examples of "twist fields" . Twist fields exist in 
a QFT model whenever there is a global internal symmetry a (a symmetry that acts 
the same way everywhere in space, and that does not change the positions of fields): 
J dxdr C[aip]{x,T) — J dxdr C[ip\{x,T). In the model with lagrangian there IS a 
symmetry under exchange of the copies. The twist fields defined by Ijl7p . which have 
been called branch-point twist fields jl7| , are twist fields associated to the two opposite 
cyclic permutation symmetries i i—^ i + 1 and i + 1 i {i = 1, . . . ,n, rt + 1 = 1). We 
can denote them simply by 7^ and 7^, respectively 

Tn = To- , a : i 1-^ i + 1 modn , (19) 

7n = , : i -I- 1 1-^- z mod n . (20) 

Notice that 7^ can be identified with 7L„. Thus for the n-sheeted Riemann surface 
along the set A made of N disjoint intervals we have 

^K„,„ « {%{ui,0)fn{vi,0) ■ ■ ■ Tn{uN ,0)fn{vN ,0)) cl.^) c ■ (21) 

This can be seen by observing that for x € [uj,Vj], consecutive copies are connected 
through T = due to the presence of Tn{vj,0), whereas for xinB, copies are connected 
to themselves through t ~ because the conditions arising from the definition of 
Tn{uj,0) and Tn{vj,0) cancel each other. More generally, the identification holds for 
correlation functions in the model C onTZ 

, , (r„(wi,o)r„(wi,o)---r„(MAr,o)r„(?;w,o)o,(x,r))£(„) c 

(0(x, r; sheet i))c n„ n — = ' i (22) 

where is the field in the model coming from the i*^ copy of £, and the ratio 
properly takes into account all the proportionality constants. 

It is often useful to introduce the linear combinations of the basic fields 

n 

^fc = ^ e^-^^-'.p, , fc = 0,l,...,n-l, (23) 
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that get multiplied by e^'^*'^/" on going around the twist operator, i.e. they diagonahze 
the twist 

In^Pk = e ' </3fc , i„<Pfe = e ' i,5fe . (24) 

Notice that when the basic fields ipj are real then (p^, = (pn-k- When the different 
values of k decouple, the total partition function is a product of the partition functions 
for each k. Thus also the twist fields can be written as products of fields acting only 

on (pk 

n—l n—l 
fc=0 fc=0 

with Tk^n'fk' = <^fc' if k ^ k' and Tk^n'fk = e^^'^^^'^ipk- In terms of these fields, the 
partition function on the initial n-sheeted surface is 

n-1 

Zn„,N = ]^(Tfe^„(Mi,0)Tfc,„(?;i,0) . . .7;i(wAr,0)T„(wAr,0))£(„) c • (26) 

fe=0 

This way of proceeding is very useful for free theories when the various fc-modes 
decouple leading to Eq. (p6| . However, as we shall see soon, this is not straightforward 
for the problems we are considering here, because the compactification condition 
introduces a non trivial coupling between the /c-modes. 



3. Compactified boson on a Riemann surface 

We consider a complex bosonic free field with Euclidean Lagrangian density 

C=j^\Vv\\ (27) 
The replicated theory is then {ipj — (pj^i + 'i^j,2) 

n n— 1 

j=l fc=0 

with partition function 

n—l \ I n—l 



= il / Wk][d(pl]e^v\-^ [ {d:,^kd,^l + d,rkd,^k)d^z\ , (28) 

where C stands for the restriction conditions in (jl8p and Ck are the corresponding 
conditions on the fields (pk- 

The fields ^pj are free, but with each component compactified on a circle. Since 
the field is complex, the target space of £ is a torus with radii Ri and i?2 ■ Encircling 
the branch points can lead to a non-trivial winding that can be written as (in the case 
of a branch point at the origin) 

ipj{e^^^z,e~'^^^z) = ipj^i{z,z) + Rimj^i + iR2mj^2 , iTij^i,m.j^2 G Z, (29) 
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where Z is the set of integer numbers. In the following we will only consider equal 
compactification radii R = Ri = i?2 • For the fields ipk these conditions read 

n n 

^k{e^^'z,e-^^'z) = e^^'^^k{z,z)+Rj2e^"-'rn, = ^^(^^(z, ^) + i? ^ , (30) 

i = 1 i = 1 

where mj G Z + iZ and we introduced 9k = e^'^*" . Consequently, for a given fc ^ 0, 
the target space for the fields ifk is compactified on a complicated "lattice" RA^/m 
where Ak/„ is the two dimensional set of vectors given by the generic combination 
through integer coefficients of the vectors {9k, 9^, . . . , 6*^} (we recall that 9^ = 1), i.e. 

Ai |g = 9lmj ; e Z + izj . (31) 

This complicated structure of Aj,/„ defining the space where the fields ipk are 
compactified is the precise reason why the approach based on Eq. (j26p becomes 
difficult for the calculation of Z-ji^^^. Few simple remarks on the structure of Aj./„ 
are in order. For fixed rt, Afe/„ = Ai_fc/„. The vectors 9\. are not independent (e.g. 
YJ]=l^i = 0). In the following we will mainly use the results by Dixon et al. |I3j . 
but we have to stress here that there are a series of papers from late eighties about 
conformal field theories on orbifold (e.g. [TH [181 [191 [20] ) that are very useful for the 
problem at hands and even for more complicated cases. The strategy of Ref. [13] 
to calculate these partition functions consists in splitting the field Cpk in a classical 
and a quantum part Cpk — (^^' + (^^". The evaluation of the Gaussian functional 
integral naturally divides into a sum over all classical solutions, times exponential of 
the quantum effective action. The classical solution takes into account the non-trivial 
topology of the target space as 

^^.'(e2--z, e'^^^H) = 9i(pt{z, z) + « , (32) 
where v € RAk/m while the quantum fluctuations are transparent to it 

^Tie'""''^. e-2-^-z-) = 9i^ri^^ ■ (33) 

3.1. The two-point function and entanglement of a single interval 

The two-point function of the twist fields for a complex field cp is il3j (see also 
[Appendix A[ ) 

(Tk^n{u)fk,n{v)) CX \ . , (34) 

\u — v\ ^^'^ 

where the dimensions of the twist fields read 

1 k / k\ 

Using Eq. (j26p . the partition function on 7^n,i is 

n—l 71 — 1 ^ 

^K„,i = n = n^'^fc."^")'^^'"^'")) " I _ ^4^„ , (36) 

k=0 k=0 ' ' 



with 



k=0 



Entanglement entropy of two disjoint intervals in CFT 



9 



and we defined the normalisation constant c„ according to [7]- This is in agreement 
with the direct calculation [7| for central charge c = 2 because we are dealing with a 
complex field. For a real field the correlation function is the square root of the previous 
result, leading to an exponent a;„ that is the half of above, in agreement with a c = 1 
theory. 

In the case of the two-point function it is then very easy to use Eq. ([26|l to obtain 
Ztj^ J. This is not the same for Z-ji^^^ and we will adopt a shghtly different strategy. 



4. The four-point function of twist fields 

Using Eq. (j26p one is tempted to write the partition function on TZn.2 as 

n— 1 n— 1 

,n{vi,0)Tk^n{u2,0)%,n{v2,0)) . (38) 

fc=0 fc=0 

This choice is consistent with the requirement that there should be no monodromy on 
going around both ui and vi or both U2 and V2- 

Using global conformal invariance the four-point function of twist fields can be 
generally written as 

7 f X f \U1 - U2\\V1 ~ V2\ Y^"^" ^ I . .„o^ 

^fe,n(Wl,fl,M2,f2)0C -j r Qk,n{x) , (39) 

V |M1 - Vl\\U2 - V2\\Ul -V2\\U2-Vl\J 

where x is the four-point ratio given in Eq. In the case of interest here x is real, 
but in the most general case can take any complex value (see also [H]). From Eq. 
((55)) Zqi^ ^ reduces to Eq. ([5]), after we normahse such that ^n{^) = 1 (for a; ^ 
Z-ji^ 2 is the product of two two-point functions previously calculated and normalised 
with c„). 

However, because of the compactification conditions of the field (^fc Eq. ([30| . it 
is not easy to write the classical contributions to Zk,n and to avoid multiple-counting 
of the various classical solutions when summing over them. For this reason, we will 
adopt a mixed strategy. We will use Eq. ([55]) for the quantum part of Z-ji^ .^ that is 
transparent to the compactification conditions. This will allow to re-use the results 
of Ref. |13| without modification. For the classical contribution instead, we will sum 
over all the possible configurations of the fields Lpj that have easier compactification 
condition than ipk. 



The quantum part 

The quantum part of the four point correlation (j39p is independent from the 
compactification of the target space and it is responsible for the full scaling factor. It 
has been calculated in Ref. [TSj and we report its derivation in the [Appendix A| The 
final result is Eq. (|A.30|1 that for real x becomes 



\U1 - U2\\V1 - V2\ \4Afc/„ ^ 



(40) 



\ui - Vi\\u2 - V2\\ui ~V2\\U2-Vi\j h/n{x) 

where 

Ik/n{x) = 2Fk/Jx)Fk/n{l -X)^ 2(3k/n[Fk/n{x)]^ , (41) 

where Fy{x) and Pk/n{x) are given in Eq. (fTU)) . Compared to Ref. [T3] and to the 
appendix we have been stressing all the dependence on k/n of the various functions. 
Note that Eq. (|4ip is manifestly invariant for x 1 — x. From this expression, the 
contribution of the quantum fluctuations to Zn^ ^ is readily obtained from Eq. ((38)) . 
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4-2. The classical part 



The value of the action on the classical solution of the equation of motion contributing 
to Z^' has been derived in Ref. [T3] for a general orbifolded theory. This derivation 
is reported in full details in the [Appendix A[ Here, for simplicity in the calculation 
and for physical reasons, we specialize to the case of x real. The action for a given 
classical configuration can be read from Eq. (|A.39P in the appendix (fixing a^/n = 
and using our normalisation for the action): 

2g7rsin (7rfc/n) r, ,t t , ,on , ^ 

^c, = ^ ] ' (42) 



where $1,^2 G R-^k/n are generic vectors of the target space lattice RAk/n- At this 
point we could calculate .^fe^„ by summing over the vectors in the target space i?Afc/„ 
given by Eq. (^1]) . This computation would be very difficult. But not only, if in fact 
we would have been able to make this sum, the total partition function Ztz„ 2 would 
not be given by Eq. ([55]) because some classical solutions would be counted more than 
once. 

For all these reasons, we prefer to calculate the total partition function Z-ji^ ^ as 
the sum over all the classical configurations independently from the value of k. Thus, 
for a complex field compactified in both directions with the same radius, we have 

Z^..- E 'liwln^ (43) 

meZ2" k=0 

and so in the total J-n{x) the scaling factor simplifies to give 



E 



const 



feOb Pk/n[Fk/n{x)]' 



exp ■ 



2.g7rsin U^) 



\^l\'Pk/r. 



I6P 



Pk/r. 



(44) 



Notice that we keep J2m '-'^ Yik because, for each fc, all the components of the 
vector TO £ Z^" (we are dealing with the complex (pj's; for real (pj's m G Z") are 
involved in the condition pop . The quantum part does not depends on m, and so we 
can factor it out as anticipated above: 



Given 

we have 



n 



const 



mgZ2" k=0 



n '^^p 



-2g7rsin (tt^) 



n-1 

i?E^i(™u 

1=0 



im 



1,2 ) ' 



1,2, 



(45) 



(46) 



E 

.«=1,2 



m'f]^ cos 



27r-(r 

n 



I / ip) (p) 

+ (™r,l'"sT2 



(P) 

m ], I TO, 



) sin 



ZTT— (r — 5) 



E[^ 

Li3=l,2 



,(p) 



(47) 
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where m'f'' G Z" and 



k 

27r-(r - 
n 



27:— (r — s) 
n 



(48) 



Notice that C k is invariant for k <-> n—k, while Sk_ changes sign. For this reason, when 
summing over k in the total partition function at fixed r, s, the two terms with Sk_ 
and S n-k cancel out. This fundamental simplification does not happen if we consider 

(rt) 

only the partition sum at fixed k. The remaining two sums over rriq € Z" factorize 
and we have 



E n 



R sin TT 



.meZ" fc = 



k 



exp I i TT 



• • TO + m' • ■ 



l3k/n m*- -C^ -m, 

}1 



H 2 



where the matrices £7 and are 



. n-l 



sm TT 



. n-l 

-E 



(3 k cos 



/3ic 



27r— (r — s) 
n 



^ k , 

2n—(r — s) 



(49) 



(50) 



(51) 



and the indices r and s run over 1, . . . ,n. We remark that the term corresponding to 
fc = is zero, thus we can deal with X]fe=i these definitions. All elements of these 
matrices have vanishing real part. 

Given a G x G symmetric complex matrix F with positive imaginary part, the 
Riemann-Siegel theta function is defined as in Eq. (fTTj) . However, H. and fl have one 
common eigenvector with vanishing eigenvalue that is (1, 1, . . . , 1). The sum in Eq. 
(|49p is then divergent and cannot be written as 0(O|il)0(O|ri). After diagonalizing 
these matrices, we will see that is easy to adsorb this divergence in the normalisation 
factor. ^ 

The eigenvalues of and SI are 



2gR sin I tt— )j/3g/„, u}q = 2gR sin I n— 



nj Pq/n ' 



9 = 1,. 



. ,n, 



(52) 



For q — n (or equivalently q = 0) the eigenvalue is vanishing, and then the imaginary 
parts of fl and ft are not positive definite. Notice also that ujq = tOn-q and ujq = Qn-q- 
The matrices Q. and have a common eigenbasis whose normalised eigenvectors 
can be written as 

Jl-Ki — r 



{yq)r = 



g,r = 1, , 



Thus, the n x n complex matrix U whose elements are Ur 
simultaneously diagonalizes Q and SI, i.e. 



unw = 



UJl 










/ OJi 















unu'' = 















) 






\ 





/ 



(53) 

{yr)s is unitary and it 
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In order to extract the divergence from (j49p . we introduce a regulator by setting ie 
(0 < e <C 1) instead of for the last eigenvalue i.e. 



(54) 





/ UJi 







] 










\ 












































[ 





ie 


) 




v 








^e J 



We introduce U as the restriction of U to the first n — 1 eigenvectors (i.e. we have 
dropped the eigenvector generating the kernel) Uqr = Uqr iov q,r — 1, . . . , n — 1. We 
remark that U is not unitary. We write the vector of integer numbers m as to = M+m, 
with M belonging to the kernel (i.e. proportional to (1, . . . , 1)) and rh to the space 
orthogonal to it. Using the orthogonality of M and rh, we have 



TO* • • TO = {Um) ■ fid ■ (Um) = [ lim^ {UMf ■ fl^^e ■ (UM)] + {Umf ■ ■ {Urn) 

where we explicitly used that in the space orthogonal to the kernel the product gives 
a finite result (and then the sum will be finite). We have the same relation for f2. We 
can re-organize the sum in Eq. (|49p . summing before on the numbers spanned by rh 
(that give a finite function in an n — 1 dimensional space) and after the diverging 
sum over the kernel. We finally find 

n 2 



Zr] — 



lim — 



where we defined 



r]T = 



I 



\ 



^1 



e(o|ryr) e(o|?7r) 



r]Y = C/t 



(55) 



(56) 



^n-l 



V 



which are symmetric and have positive imaginary parts and therefore provide well 
defined Riemann-Siegel theta^ functions. _ 

The matrices r]T and r]T are D, and fl respectively with the last line and row 
dropped, i.e. 



2i 



2i 



ji-i 



k = l 
n-1 



^ sin ( TT- j (3k/n COS 



sin ( TT 



fe = l 



ZTT— (r — s) 

n 



27r— (r — s) 
n 



fc=i ^ ^ 

n— 1 ^ 

sin ( TT 



if3. 



k/r, 



fe = l 



„277ii(r-s) 



„27ri|.(r-s) 



where r, s = 1, . . . , rt — 1. We introduced the matrices T and F, in such a way that 
they do not depend on R and we defined 

V^9R\ (57) 

that is exactly the same as in Ref. [5], while the normalisation of R is different. 

Thus, taking pS)) and (|55p . adsorbing the divergence for e ^ into the constant, 
we have 



const- 



[e(o|7?r) 9(o|77r)]^ 

lTk=\Pk/n[Fk/nix)]' 



const 



[e(o|77r)9(o|r?r)]^ 

nfc=\ Fk/n{x)Fu/n{^ - X) 



(58) 
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This expression is manifestly invariant for a; — > 1 — a; because under this transformation 
Pk/n ^ l//?fc/n and so r ^ f . 

In order to fix properly the normalisation constant (given by J-"„(0) = 1) and to 
show explicitly the invariance under ry — > 1 /?], it is worth to manipulate this expression. 
Using Poisson resummation formula, in [Appendix B| we show the following identity 



e(o|ryr) = T]- 



Y[f3k/n) e(o|r/77) 



(59) 



We finally have 



const- 



e(o|77r) e(o|r/r?) 



(60) 



Uk=l(^k/n [Fk/n{x)] 

where we fixed the constant by requiring that for J-n{0) = 1 ( we used that for a; ^ 0, 
-Pfc/n(0) = 1, and /3fc/„ +oo; furthermore 8(0|r) goes to 1 for /3^./„ ^ cx)). We have 
thus written J-n in a way that is manifestly symmetric for the exchange ry <-!■ l/ry. 
A last manipulation can be done by using 



[0(o|r)p= n ^k/nix), 



(61) 



fe = 1 



proved in [Appendix C| This finally leads to the square of Eq. ([5]), in fact this equation 
is valid for a complex field compactified in both directions with the same radius. The 
real field corresponds to the square root of the previous result. This final manipulation 
allows to manifestly show that J^nix)\r]=i = 1 identically. 

It is also worth to mention that we wrote the partition function over a n-sheeted 
Riemann surface in terms of a Riemann-Siegel theta function defined with a matrix 
of dimension n — 1 that is exactly the genus of the covering surface. 



4-3. Special cases 

4-3.1. n ^ 2. In this case the matrix F is just 1 by 1, and so 8(0|r) is a standard 
Jacobi 03 function. Thus Eq. (l60|) (ti/2 ~ iPi/2) 

%{ti/2vW3{ti/2/v)'''^ 



^^(Tl/2) 



(62) 



where we used that ^1/2(2;) = G\(t-yI2)- This is exactly the square of the result in Ref. 
[9] as it must be. 



(63) 



4.3.2. n — 3. First we observe that there is only one t because ri/3 = T2/3. We have 

r = liZi r 2 = J_f ^ -1 \ / 7/3 -1 

V3 I -1 2 ; V2 I 1 1 A 7 J V2 I 1 1 

where 7 = \/3ti/3. Since these matrices are written in terms of integers only, this 
allows to write 0(O|r7r) as a sum of ^3 and 62 and finally using the duplication formulas 
to prove the following identity 



92(^7)^^2(f)' + 03(^7)^e3(f)'+^4(r77)^04(f^' 



(64) 
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The same is clearly true for 9(0|r/r/), obtaining 
1 

.77/ Vdry/ \6r]/ yrjJ \6rj 



2 



4.3.3. n = A. Using tricks similar to the case n = 3, it is possible to write 0(O|r) 
for n = 4 as a sum of 6 terms that are products of three Oi (with different r's now, 
because in general T2/4 ^ T1/4 = T3/4). This expression is however very long and not 
illuminating. 

4.4. Decompactification regime 

For fixed x, in the limit of large 77 we have 0(O|7/r) = 1 + . . . and 
where . . . denotes vanishing terms as ry 00. Thus 

J-n{x) — , (66) 

nfe=l Fk/nix)Fk/n{l - X) 

recovering the correct quantum result Eq. (|A.30[) with the proper 77 dependent 
normalisation (we used det(-ir) = JJ^Zl Pk/n and (3k/n = Fk/nC^ - x)/Fk/n{x)). 

Note that, using the symmetry 77 <-> 1/77, the same formula gives also the answer 
for 77 ^ 0. 

4.5. Small X regime 

Another case when the final result can be written in terms of simple functions is the 
expansion for small x, corresponding to two far distant intervals (and by symmetry 
a; — > 1 — a; also for x close to 1, corresponding to close intervals). 
For X —> the expansion of /3fe/„ is 



where "fE is the Euler gamma, ip{z) = r'(z)/r(z) is the Polygamma function (also 
know as digamma function) and pi{z) is a polynomial of degree 21, whose explicit 
expression is not needed. Plugging this expansion in the Riemann-Siegel theta function 
we obtain 

e(o|77r) = 1 + (67) 

x'li T>l=l ''H'^iT Cfc/„- mg,,f j2Zl M^'iT f"/^ c'fc/„- ™ (-1 _^ o{x)) , 

meZ"-i\{0} 

where the matrix Cj./„ is defined in (|48p . In this expansion the leading term is provided 
by those vectors m S Z"~^ \ {0} which minimize the expression 
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This expression is obviously minimized by all vectors of the form m\ = 
(0, . . . , 0, ±1, . . . , ±1, 0, . . . , 0), with I contiguous ±l's . At fixed I = 1, . . . , n — 1, 
there are 2{n — I) of such vectors for which the expression ([55)1 is 1. In order to 
evaluate the coefficient in front of the leading term, we observe that 

Using the following integral representation for the digamma function 

-* _ p-yt 



V'(2/) + 7i=;- / n —"^^^ ^^^^ 

exchanging the order of sum and integral and then performing the latter, after simple 
algebra we find 



2 / fc^ ^ 



7r- I /fc/„ • Ck/n ■ mL± = - log 



2n sm TT— 



2 



(71) 



— > sm 
n ^ 

fc = 1 

Now, assuming 77 ^ 1, in the ratio of Riemann-Siegel theta functions occurring in 
J>i(.t) the leading term is given by the minimum between 77 and 1/77, that we indicate 
as a = min(77, 1/77)) Therefore we get the small x behavior of the scaling function 

n — 1 ^/ ,^ n— 1 , 

^„(..) - 1 + X" E + ■ • ■ - ^ + Ki^) ° E 7-7^ + ■ ■ ■ ^ (^2) 

(=1 [2"sm(7r^)] / [sm(7r:^)] 

where the dots denote higher order terms in x. Note that for n — 2 this results reduces 
to Ti(x) = 1 + 2(a:/16)"""('''^/''^ + . . ., as already found in [S]. We stress, once more, 
that this expansion is valid only for 77 ^ 1, in fact for 77 = 1 the denominator in Tn{x) 
(that is of order 0(a;)) cancels exactly the numerator. 



J1..6. Different compactification radii 

In this manuscript we mainly condisider the case of a complex boson compactificd in 
both directions with the same radius, because it has more physical applications and 
also to lighten the notation. However it is straightforward to generalize to the case 
with different compactification radii (at least for real four-point ratio x). The only 
change is that the target space is given by the product of two circles with different 
radii i?i 7^ i?2- Eq. now becomes 

n-1 

= Y.(^i{Ri^tl+^R2^ti) ■ (73) 

(=0 

The only (minimal) changes compared to the case of equal radii are in the computation 
in Eqs. ([T7|) . Repeating the straighiforward algebra we have 

n-1 
r,s — 

n-1 

= E [(^1),, (Rl^rWl^l + Rlmi^mi^}) ~ R1R2 (S^)^^ {m^^W^l - m^^lm^^^) 

r,s — 

= [ Rj m[''^ * • • m[P'> + Rl m'f'' * • Ci • m^"' + R1R2 tr (a 5^) " . 
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The important point is that the term in R1R2 is still vanishing and so the classical 
part of the action corresponds to the one of two independent real bosons. Note that 
this property is not completely trivial because the compactification conditions are not 
on the fields (fk- The divergences due to = are removed as before, arriving to the 
final expression 



e(o|77ir)e(o|r/7?i) 



e(o|r)2 



e(o|?72r)e(o|r/772) 



(75) 



where rjj = gRj, the matrix T is defined in Eq. 
chosen such that ^n(O) = 1. 



e(o|r)2 

^ and the normalization has been 



5. The analytic continuation and the entanglement entropy 

In order to obtain the entanglement entropy we should be able to analytically continue 
Eq. ([5]) to general complex value of n and, only after, take the derivative for n ^ 1. 
The definition of the Riemann-Siegel theta function makes this program hard, because 
it explictly involves summations involving a matrix of dimensions (n— l)x(n— 1) 
that are not obviously continued to real values of n. One should find a different 
representation of the same function which is manipulable and we have been unable to 
doit. However, it is possible to analytically continue the denominator of Eq. ([5]). This 
is not only a first step towards the full analytic continuation, but it provides all the 
answer in the decompactification regime (Eq. [BB]) , allowing to give precise predictions 
for <C 1 and 77 1. 

The logarithm of the denominator in Eq. (|60p is a sum over fc = 1 . . . n — 1, which 
involves analytic functions of fc. In this case it is possible to use the residue theorem 
to write the sum as an integral in the complex plane of logF^/ni^) times a function 
that has poles only for integer values of k and over a contour that encircles all of them. 
A useful representation of this is 

Dn{x) = y^logFj./„(x) = — 7rcot(7rz)logF2/„(a;) , (76) 



k=l 



where C can be chosen as the rectangular contour in — iL, n + iL, iL, —iL) (because 
logFo(a;) — logFi(x) = and cot(z) has no poles for Imz 7^ 0). We can now change 
variable z/n — *■ z to obtain 



Dn{x) = n 



f dz 

/ — cot(7rzn)logF,(x), (77) 

and C is the re-scaled rectangle (1 — iL, 1 + iL, iL, —iL). This simple formula provides 
the desired analytic continuation. Notice that all the poles of the integrand in the strip 
< Rez < 1 are on the real axis. Then, if the argument of the integral would decay 
quickly enough for Imz —^ ±00, we could send L —^ 00, ignoring the contribution of 
the horizontal pieces and remaining only with the vertical ones. This is unfortunately 
not the case, because the integrand is increasing when L 00. 

We can however take the derivative wrt n in the contour integral. This leads to 

dDnix) 



D[{x) 



dn 



n=l 



I ^-^logF,(x). (78) 
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Figure 2. Mutual information I^^.j^^ for tlie XXZ model. All numerical data are 

extracted from Ref. [9]. Left: for x = 1/2 (top curve) and x = 1/4 (bottom 

curve) as function of r]. The continuous curve is the decompactification result Eq. 
I|82| l. Right: at fixed rj as function of x. The top curve corresponds to 

7] = 0.295 small enough to agree for all considered x with Eq. H82| l. The other 
curves correspond to higher values of rj, when the small r] approximation looses 
validity. 



In this integral, the horizontal contribution is vanishing for L — > c» and so only 
the vertical ones are left. Because of the periodicity of the integrand these two 
contributions are equal and so 

^-^logF.(x). (79) 

-joo * sm -KZ 

Such integral is easily evaluated numerically for any x. For x = 1/2 it is possible 
to get an analytic result with a different method (see [Appendix D[ ) that agrees with 
the value calculated numerically. To cross check our results in the appendix we also 
provide the analytic continuation as perturbation series in x. Despite the asymptotic 
character of the perturbative expansion, it provides a very good approximation for all 
X < 1/2. 

5.1. The entanglement entropy in the decompactification regime 

From Eq. (|66p and using the result above for the analytic continuation, we have that 
the entanglement entropy for a real boson in the decompactification regime is 

5.(,«l,-S?4l„,-2<i>±£i(i^, ,80) 

where is the result in Ref. 7\. The same result obviously holds for 77 3> 1 with 
the replacement 77 ri^^ . 

This should be compared with the numerical results for the XXZ model by 
Furukawa et al. [9], where the entanglement of the XXZ chain for generic values of 
the anisotropy A and magnetic field (always in the gapless phase) has been calculated 
by direct diagonalization for systems up to 30 spins. In the absence of the magnetic 
field, rj is related to the anisotropy by 77 = 1 — (arccos A) /tt, while for non-zero a 
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Figure 3. In in Eq. l|85|l as function of rj for x constant x = 1/2 (left panel) 
and X = 1/4 (right panel). In each plot the three curves corresponds to different 
values of n = 2, 3, 4 (from bottom to top). 



closed formula for 77 does not exist and must be calculated numerically as explained 
in [9] . The main results have been written in terms of the Rcnyi mutual informations 



r(") _ c(") I c(") 
'Ai:A2 — "^Ai "t" ''As 



S 



in) 

A1UA2 



1 



clog(l - x) + 



1 



(81) 



6n n — 1 

where Ai and A2 are the two intervals composing A = Ai U A2. For n = 1 the 

data (reported also in Fig. [2|), when plotted in terms of x collapse on a single curve, 

(2) 

confirming the validity of the scaling form. In Ref. 9 , only has been compared 

with the available CFT prediction Eq. ([7]). However, for n = 2, the collapse of the 
data is worse than the one for I'^^.j^^ because of the strong oscillating corrections to 
the scaling of the Renyi entropies (analogous to the ones observed for a single interval 
[23l[22]). However, the agreement was rather satisfactory, considering the small system 
sizes and the oscillations. 

Here we are in position to offer a first prediction for the von-Neumann mutual 



(1) 

Ai:A2 



The numerical data from Ref. [9] are reported in Fig. [2] The 



information / 
prediction in the decompactification regime is 



(1) 

Al-.A; 



_l.,D[{x)+D[{l-x) 
2 ' 2 



(82) 



where again l')^^'^^ is the result of Ref. Jj . This prediction, for various values of 77 and 
X are reported in Fig. [21 In the left panel it has been reported the mutual information 
as function of rj at fixed x = 1/2, 1/4. On the scale of this plot, the decompactification 
prediction reproduces the numerical data up to 77 ~ 0.4 and for larger values clearly 
deviates. In the right panel of Fig. [2] it is reported the x dependence of the mutual 
information at fixed rj. Again for the smallest value 77 = 0.295 the decompactification 
approximation is valid for all x. Increasing 77 this is no longer true, but it works better 
for X ~ 1/2, while it quickly deteriorates moving from the central point. 



5. 2. The entanglement entropy for small x 

We can obtain the analytic continuation also in the small x regime starting from the 
results in Sec. 14.51 We need to analytically continue Tn{x) in Eq. ([72)) valid for 77 ^ 1. 
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Figure 4. /„ in Eq. II85I I as function of x at constant rj = 0.3 (left panel) and 
r] = 0.5 (right panel). In each plot the three curves correspond to different values 
of n = 2, 3, 4 (from bottom to top). 



Calling a — imn{ri, l/rj), we can write J-n{x) as 

n—l , , 

.F„(x) = l+2n(— ) + with Pn = Y. ,. ,i2. -(83) 

Within this definition we have that, after analytically continuing, the contribution to 
the entanglement entropy is (we recall that Pi — 0) 

Sa-S^ ^ -J"i{x) = 2i-2«a."p| + . . . . (84) 

The most important consequence is that for small x the entanglement entropy has 
a power law behavior in x with an exponent that is always a for any rj ^ 1. 
The multiplicative coefficient of this power law is exactly calculated by analytically 
continuing P„. This derivation is however rather cumbersome and it is reported in 
[Appendix D.3 The dependence of P[ from rj can be read in the plot in Fig. IDll 



5. 3. Results for integer n 

In this section we report some explicit results for integer n. These are shown in Figs. 
[31 [H and [HI where we always plot the quantity 

In = log J^„(x) , (85) 

n — l 

that contributes directly to the Renyi entropy and mutual information. After analytic 
continuation, this quantity has also a smooth limit to n = 1, and so it is ideal to show 
some general properties. 

In Fig. [21 following Ref. [S] (also left panel in Fig. ^ , we plotted /„ as function of 
rj with X kept constant to x = 1/2 (left panel) and x (right panel). In each plot 

the three curves correspond to different values of rt = 2, 3, 4. It is evident that on the 
scale of the plot, the differences between various n are tiny. This means (if nothing 
really strange happens in the analytical continuation) that also the equivalent plot for 
Ii (contributing to the entanglement entropy) will be qualitatively and quantitatively 
similar. In particular, to appreciate the differences between various values of n on these 
kinds of plots, the numerical results must be extremely precise. In fact, the results 
in Ref. f9 (reported also in the left panel of Fig. [2|) are practically indistinguishable 
from those in Fig. [3l 
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Figure 5. In in 
Eq. I I85I I as func- 
tion of n at con- 
stant T] and X. 
The points at n 
integer are real 
data, the lines are 
only guide for the 
eyes, except the 
top one that is the 
analytic continua- 
tion for small r). 



A more effective way to show the differences between various n is to plot /„ as 
a function of x keeping 77 constant, as in Fig. [H Because of the reduced scale of 
variation of /„ with x, the differences between the various n are evident and guessing 
quantitative features of the analogous plot for the analytic continuation at n = 1 is not 
recommended. When n increases the corresponding value of /„ increases, in qualitative 
agreement with the numerical results of Ref. |9] (see Fig. 3 there). As already stated, a 
direct comparison of the numerical and analytically results is quantitatively impossible 
because of the strong oscillating corrections to the scaling in the XXZ model. 

Finally in Fig. [5] we keep constant both 77 and x and show the dependence of In 
on n. On the scale of the plot, the n-dependence is small. However, before making 
wrong considerations about a smooth behaviour in taking the limit for n — *■ 1 in 
the analytic continuation, it is worth to have a look at the exact results we have 
for small t]. We then consider rj — 0.2 and x = 1/2 (top curve in Fig. [5|). The 
points are the exact results for n integer from Eq. ([8|), while the continuous line is 
the analytic continuation given by Eqs. (jP.Sp and (|66p. For such small value of 7], 
the decompactification formula is indistinguishable from the exact data. The analytic 
continuation displays a pronounced binding of the curve for 1 < n < 2, that makes 
difficult any naive extrapolation from the data with n > 2. 

6. Conclusions 

We considered the entanglement of two disjoint intervals A — [ui,vi] U [^2,^2] in 
the ground-state of the CFT of the Luttinger liquid (free compactificd boson). We 
presented a complete analysis for Trp^ for any integer n, leading to the result we 
anticipated in the introduction Eq. ([4]). We have not yet been able to analytically 
continue the numerator of this expression to obtain the entanglement entropy. 
However, in the decompactification regime when the exponent rj is very small or very 
large, we calculated explicitly the entanglement entropy and the associated mutual 
information leading to Eq. (|82p . Our predictions agree well with the numerical 
computations of Furukawa et al. [3] for the XXZ model that is described by this 
CFT. 

The error in Ref. [7] also carries over to the case of semi-infinite systems, where 
now half of the points are the images ones. Also in this case a still unknown function 
J^n{x) corrects the result. However, this error does not impact on the later application 
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by two of us of these methods to studying the time-dependence of the entanglement 
entropy following a quantum quench [32] . This is because in this case only the limiting 
behaviour, when the various cross-ratios are either small or large, is needed, and this is 
insensitive to the precise form of J-n{x). Likewise, as far as we are aware, all the other 
conclusions of Ref. [32] remain valid. Instead the function J-n{x) affects some results 
in Ref. |33j for local quantum quenches, that however are easily corrected using the 
formulas for general four-point correlation functions presented in the same paper. See 
for a more detailed discussion the review [3]. 

The analytic continuation of to general real n remains the most interesting 

open problem. However, there are several other issues that deserve to be discussed. 
Firstly it is difficult to check directly our predictions for Trp^ for integer n with 
the numerical data in Ref. [9]. This because, as observed already in [9] for n = 2, 
there are strong oscillating corrections to the scaling making the comparison hard, 
if not impossible, for the system sizes accessible by exact diagonalization. There are 
several possible way-outs to this problem. One could try to describe these corrections 
analytically, as done for the single interval [22] and adding them to the leading 
contribution, but this seems to be very difficult. Alternatively one could use different 
numerical methods to access largest system sizes. Density matrix renormalisation 
group (eventually in the recent version proposed in Ref. [15) to deal with a similar 
issue) could be effective. Furthermore, the Monte-Carlo based approach by Caraglio 
and Gliozzi [TTj appfies to the case of integer n. This could then be used to test 
our predictions for models showing smaller oscillations (even in the same universality 
class). One can also wonder whether Eq. ([4]) is enough to calculate analytically the 
full spectrum of eigenvalues of the reduced density matrix. In fact, in Ref. [S^ it has 
been shown that this can be calculated by knowing Trp^ for integer n only. 

A few comments are also in order. In Ref. [24l |25l [26l [27l [29] , the entanglement 
entropy for two disjoint intervals has been calculated for free fermionic theories, 
that after bosonization always correspond to a compactified boson with rj — 1/2 
[12] . However, it has been found that the entanglement entropy is given by Z^^^, 
apparently in contrast with the numerical calculation in Ref. [9] and what found 
here. The details of this apparent disagreement are still not completely understood, 
but they should be traced back to the different boundary conditions that result from 
constructing the reduced density matrix for spin or fermion variables. For the Ising 
model numerical computations [28] also show a good agreement with ^ . Also 
in this case, it is likely that the deviations from Z^^ ^ should be attributed to the 
choice of the variables used in constructing the reduced density matrix. (In fact, the 
calculations in the spin variables [35j show numerically and analytically that Z^ ^ 
is not correct.) Finally holographic calculations in AdS/CFT correspondence [30l [31] 
considering the classical limit in the gravity sector, also found Z^^ ^ . It would be 
interesting to understand how the correct result might arise from taking into account 
the quantum effects on the gravity side. 

We close this paper by discussing the case with N > 2 disjoint intervals. As far as 
we are aware there are no firm results in the literature. By global conformal invariance 
we have 



For J-n,N{{x}) = 1 this is the incorrect result of Ref. [7] (note a typo in the 




(86) 
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denominator), {x} stands for the collection of 2N — 3 independent ratios that can be 
built with 2N points. Some old results from CFT on orbifolds in Refs. [TU [T^ could 
be useful to calculate Tn^N{{x}) for a compactified boson. 
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Appendix A. Correlation functions of twist fields and Z„ orbifolds 

In a order to make this paper self-contained, in this appendix we review the results 
of Ref. [13] for CFT on orbifold spaces that we used. In particular we will show how 
to obtain Fqs. (|40p and (|^^ . As we have stressed in the main text twist fields exist 
in a QFT whenever there is a global internal symmetry. An orbifold is obtained by 
identifying points of the target space through a given equivalence relation, leading 
automatically to a global symmetry and so to the presence of twist fields (that is 
why Ref. |13| contains the main ingredients for our calculations). To be explicit a D 
dimensional orbifold R^/S' is obtained by identifying points of the target space R'^ 
through an equivalence relation 

A' -A, if X' = eX + v = gX, g^{9,v), (A.l) 

where is a rotation and w is a vector of R^. The set of the pairs S ~ {{9, v)} defining 
the orbifold is called space group, while the subgroup A = {g = (l,f)} C S, made by 
the translations only, is called lattice A of 5. A simple example of orbifold is obtained 
by identifying points in the target space with opposite signs X ^ —X. This gives rise 
to so-called Z2 orbifolds. A Z„ orbifold is generated by a rotation 9 of order n (i.e. 
0" — 1) and its space group S is given by the pairs {9^ ,v) with j = 0, 1, 2, . . . , 71 — 1 
and V runs over an even dimensional lattice A. Among the conjugacy classes making 
up the partition of S, we distinguish the ones of the form {{l,9^vo)} {vq is a fixed 
vector of R^), which contain only the translation elements of S and therefore describe 
the winding sectors. Instead, elements like {0^ ,v) with j = 1,2, . . . ,n — 1 and v G A 
belong to classes like 

{ {9^ , 9'-vo + (1 - 0')u) ; r e Z, u e A} , (A.2) 

which describe the twisted sectors. Notice that, for each j — 1, 2, . . . , 71 — 1, there 
are many conjugacy classes of the form {{0^ , v)} with v belonging to different subsets 
(cosets) of the lattice A. Thus, as for the twist fields associated to the Z„ orbifold, 
each of them is characterized by two indices: the index j = 1, 2, . . . , , n — 1 of the 
twisted sector of the Hilbert space and another index e labelling the conjugacy class 
within that sector. The twist-fields in our problems are exactly the same that appears 
in Z„ orbifolds, but since we do not have identifications of points the index e can be 
only zero and will be ignored in the following. 

Let us go back to our main goal, that is the calculation of the four-point function 
of twist fields. Let us consider a complex field X{z,z) defined on the worldsheet given 




Figure Al. The two closed loops Ci and C2 we considered as basis. The different 
sheets are indicated as solid vs dashed lines. 



by the Riemann sphere (this is one of the (pk in the main text). The occurrence of the 
twist field in the origin tells us how the field X is rotated and translated when it is 
carried around this point in the worldsheet, i.e. 

Xie^'^'z, e-^-^'z) = e^X{z, z) + v , (A.3) 

where u is a vector of the coset. The phase rotation 9^ of X is the monodromy of the 
field. In order to make this definition meaningful, the field X must be complex. 

When X is taken around a set of points where different twist-fields are placed, 
it is rotated and translated according to the product of the space group elements 
associated with the different twist fields. The relevant circuits to consider in order 
to completely fix the boundary conditions are the ones enclosing a collection of twist 
fields with net twist zero, namely the loops along which the field X acquires no phase. 
Such paths C are called closed loops. Two example of closed loops that we will use in 
the following are reported in Fig. lAll 

Splitting the field X =^ Xd + Xqu into its classical part X^x and its quantum 
counterpart Xqu, we impose (IA.3I) by requiring that 

Xei(e2"z, = e'X,x{z, z) + v , (A.4) 

and 

X^^ie'^'z, ^ e^X^^{z, z) , (A.5) 

namely A"qu ignores the translations in the space group. Thus, for any closed loop C, 
we have 

AcXqu ^ I dz d.Xqu + <f dz 9jXqu = , (A.6) 
Jc Jc 

and 

AcXci = (j> dz d,Xd + f dz d,Xd = V , (A.7) 
Jc Jc 

where now v is not the same of (IA.4I) . but it is a vector of A depending on the twist 
fields enclosed by C. For the twist field 7fc_„ we have 

V e {l-e)A, (A.8) 

where, compared to Ref. [13j . we fixed /^j = /e^ = 0, because we have only the trivial 
fixed point 0. The requirements (|A.6[) and (|A.7|) are the global monodromy conditions. 
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Let us take 9 = e^^*'^/" (i.e. 9k in the main text) to be the rotation by an angle 
2'Kk/n for a fixed fc e {f , 2, . . . , n — f }. In all this appendix we will consider fixed 
values of k and n, and so there is no ambiguity in denoting with T and T what in the 
text we called 7fe^„ and 7fc_„ (i.e. the twist field associated to 9 and 9~^ respectively 
and that in [13] are called cr+ and cr_). We fix j = 1 in Eqs. (|A.4p and (jA.Sp . Also 
all other quantities will loose the subscripts k and n etc. We are interested to the 
following four point function 

z = (r(zi,zi)r(z2,z2)r(^3,^3)r(z4,^4)) = j [dx][dx]e-^[^^^i , (a.9) 

where the action for the complex field X reads 

S[X, X] = ^ J {d.XdgX + d,Xd,X) <fz . (A.IO) 
Separating the classical contribution from the quantum part, we can write 

^ = ^qu^e-^^', (A.ll) 

where Sc\ = S[Xc\, Xd]- The classical part in (jA.ll[) is given by the sum over the 
possible configurations of the classical field, which are characterized by (|A.7p . 

Following [13] . let us start by considering the Green function in the presence of 
four twist-fields 

^ (-ia.xa.xr(.ot(..)T(.3)r(.,)) _ 

Imposing that for z —f w we have g{z,w;Zj) ^ {z — w)^"^ and that for z —>■ Zj we 
have g{z, w; Zj) ^ {z ~ Zj)~^/^ and g{z, w; Zj) ~ (z — Zj)^^^^''/"-'> for j odd and even 
respectively (and the opposite for w — > zj), we can write g{z, w; zj) as 

■fc (Z - Zi){z - Z3)(W - Z2){W - Z4) 



g{z, w; Zi) = ujk{z)ujn-k{'w) 



n [z - 



k\ {z - Z2){z - Z4){w - zi){w - z:i) ^ ^1 _ . 
n J [z — w)^ 



(A.13) 



where 

u^kiz) = [{z ~ z,){z - Z3)]-"/"[(^ - Z2){z - Z4)]-'^-'=/") , (A.14) 

and A{zj,Zj) is a constant (in z and w) that contains the dependence of g{z,w]Zi) 
on the antiholomorphic coordinates and must be determined by global monodromy 
conditions. 

Let us now consider the limit w — > z 

r r / A / ^-2^ I k ( ^ k\ f \ 1 1 1 ^ 

\\m[g{z,w)- {z-w) J = -- 1 



2 n \ n J \z — zi z — z^ z — Z2 z — Z3 



(A.15) 



(Z - Zl)(z - Z2)(Z - Z3)(Z - Z4) ■ 

This is exactly the expectation value of the insertion of the stress energy tensor of the 
field X in the four-point correlation function. By taking the limit z Zj for any j, 
one reads the dimensions of the twist fields T and T 

n '^k f k\ , , , 
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as anticipated in the main text. Let us also introduce the auxiliary correlation function 

h{z,W;Zj) = = B{Zj,Zj)uJn-k(zMn-k{w), 

{T{zi)T{z2)T{z3)T{z4,)) 

where the rhs above comes from the same considerations as before, but the singular 
terms in (z — w) do not occur because h{z,w; Zj) is regular in this variable. 

In order to lighten the notation it is convenient hereafter to employ the following 
conformal map 

(Zl ~ Z3){Z- Z4) 

which sends zi, Z2, z^ and Z4 into 0, x, 1 and 00 respectively, where x is the four-point 
ratio {zij = Zi — Zj) 

filial, x{l-x)^ '^''f'f'^\ (A.18) 

Z13 Z2A 2:31 

After this mapping the dependence on Zi in the varius functions reduces to the 
dependence on the ratio x and its complex conjugate x only. 

Taking the limit z — > x in Eq. (|A.15P , one finds the following differential equation 
for the quantum part of the correlation function 

a.lnZ,.(.,.) . -2A,/„ - ^) - , (A.19) 

The global monodromy conditions allow to determine A{x,x). In terms of the 
functions g and h they read 

= ® dzg{z, w) + f dzh{z, w) . (A. 20) 

Jc Jc 

Dividing by LUn-kiw) and letting w 00 this gives 

A{x, x) j) dzLJk{z) + B{x, x) j) dzujk{z) = — ^1 — — ^ j) dz{z ~ x)uJk{z) . (A. 21) 

These integrals are valid for all closed loops. We choose as a basis of these loops the 
two ones depicted in Fig. lAli that suffices to determine A{x, x). 
All these integrals are easily calculated giving [T3] 

dzwfe(z) = 27rie""''/"F(a;), 1 dzujk{z) ^ 2mF{l - x), (A.22) 

dzu)n-k{z) = 2me-'''^/'''F{x), <f> dzCon^k{z) ^ -2iTiF{l - x), (A.23) 

(A.24) 



(1 - I dz{z - x)ujk{z) = 27rie-"'=^/"a;(l - a;)^^^j^ , 
V ''^J Jci dx 



1 - - ) / dz{z - x)ujk{z) = 2TTix{l - x) '^^j^^ ^ (A.25) 



n J dx 
where we introduced 

F{x) = 2Fi{k/n,l- k/n;l;x) . (A.26) 

From the two equations (|A.21[) for Ci and C2, eliminating B{x,x) and solving for 

A{x, x) we get 

A(x,x) = xil- x) r^-^ , (A.27 

dx 
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with 

I{x,x) = F{x)F{l - x) + F{x)F{l - x) = 2P{x)\F{x)\'^ , (A.28) 
and we also introduced 

t{x) EE a{x) + I p{x) EE z ^^I'^^y (A.29) 

F{x) 

For n = 2, t{x) gives the modulus of the torus, which is the covering surface of our 
Riemann sphere with four branch points, but for higher values of n this is not true 
anymore and (|A.29p has no geometric meaning. A{x, x) gives the desired quantum 
partition function 

const 1 

\x{i^^xy^^l{^x) ■ 

When this equation is specialized to x real, it gives exactly the result anticipated in 
the main text in Eq. 

For the classical part, we only need to construct the properly normalised classical 
solutions. From the equations of motion, one easily sees that dzXd and dzX^x are 
holomorphic while dzX^i and dgXd are antiholomorphic and they can be written as 

dzXci{z) = auJk{z), dzXci{z) = 6cj„_fc(z) , 

- _ ~ _ _ (A. 31) 

dzXci{z) = auJn-k{z), dgXciiz) = buJk{z) , 

where ujk{z) is given in Eq. (lAlIl) . 

The complex constants a, a, b and b are fixed through the global monodromy 
conditions (|A.7[) for the closed loops Ci and €2- (Wc do not adopt here the notation 
of [13] . where d = d and b = b, because we find it misleading.) In order to write them, 
one constructs two classical solutions Xd.i and Xc\^2 having the following simple global 
monodromy conditions 

Ac,Xei,j = Ac.Xij = 27r5„, i,j = 1,2, (A.32) 

and finds out the corresponding complex constants at, di, bi and 6^, which read 

ai = -e^""ai = -ze"" ;, _/ , 03 = 03 = -i— ^ , (A.33) 

I(X,X) i(x,x) 

9^,fcr „,kF(l — x) , r F(x) 

6i = -e^"-6i = -ie"»- _' b2 = b2 = . {AM) 

I(X,X) 1[X,X) 

Then, from (|A.7p for Ci and C2 and (jA.Sp . one gets the complex coefficients for Xc\ to 
use in (|A.31|) . which read 

a — aivi + 02^2 , b ~ bivi + 62^2 , (A. 35) 

d = divi + d2V2 , b — bivi + 62^2 , (A. 36) 

with the vectors characterizing the global monodromy condition of Xci,i and Xci,2 
given by 

vi,2e{l-9)A. (A.37) 

The last step we need for our purposes is the expression for Sd- By employing the 
following integral 



\ujk\^d^z= / — = — , , , J (x,x), (A.38) 

' ' ' ■,\^i\z-x\'i'-^)\z-l\^i sin(^fc/n) ^ ' ^' ^ ^ 
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one finds that 

where 7 = _je~''^'^/" and we introduced the vectors (j = 1, 2) independent of k/n 
that are generic vectors of the target space lattice A, having therefore Vj = (1 — 0)(,j. 

When we specialize to the problem with real x G (0, 1), some simplifications occur 
in the formulas, r is purely imaginary (i.e. a = 0) and I{x,x) — 2F(x)F(l — x). Then 
Eq. (|A.39P reduces to Eq. (|42l) in the main text, where we restored all the k and n 
dependence in each quantity. 

Appendix B. A transformation formula for the Riemann-Siegel theta 
function 

The Riemann-Siegel theta function O(0|r) for any symmetric complex matrix T with 
positive imaginary part can be re-written as 

Q{0\T) ^ J2 e""'-'^ " = /d^s J2 Sg{s-m) e"''-'^-' . (B.l) 

mGZs meZs 

Now we employ the following identity (which is a special case of the Poisson 
resummation formula) 

Plugging (jB.2p into (jB.ip and inverting J d^s and X^jezs (|B.1[) . we get a g 
dimensional gaussian integral, which gives 

Now, by applying this formula for T = AT = -A\T/X)-^A, where A = 2U''HU 
with H = diag(. . . , sin(7rfc/n), . . .), we get 

e(0|7yf) = ^det(r/(i,7))e(0|A-i(r/ry)A-i), (B.4) 

where we have used that A*' = A and dot A — 1 (here A and F are the matrix defined 
in the main text). We employ the following identity 

e{0\A-\T/ri)A-^) = e(0|r/77) , (B.5) 

which is numerically true, but we give here without a proof. Putting everything 
together we have 

1 /"-^ V^' 

©(oi^yr) = ^ [] e(o|r/r,) , (b.6) 

^ ' \k=i I 

which is exactly Eq. (|59p we wanted to prove. 
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Appendix C. A Thomae-type formula for singular Z„ curves 

In this appendix we employ a Tliomae-type formula obtained in j36j for singular Z„ 
curves to prove equation ([6T|l . In the following singular Z„ curves d,™ are 
considered 

m m 

w" = p{z) q{zr-' , p{z) - n - , q{z) = H " ^^j) ' C^-^) 

i = i = 1 

They have singularities at the points P2 = (z2,0), . . ., P2m — (•Z2m,0), while 
Pi = (zi,0), . . ., P2m+i = (2^2^+1,0), P2m+2 = Poo = (oo,oo) are the branch points. 
These curves are n sheeted coverings of the complex plane and they define Riemann 
surfaces of genus {n—l)m. The cycles ak and /S^ (A; = 1, . . . , {n—l)m) (generahzations 
of Ci and C2 in Fig. lAll for the many intervals situation, see Ref. [31] for a pictorial 
representation) provide the basis of the closed loops. We also introduce 

duj+(k-i)ni = 1 dz, j = l,...,m, fc = l,...,n-l, (C.2) 

A* 

which are a basis of the canonical holomorphic differentials. The (n — l)m x (n — \)m 
matrices A of the a-periods and B of the /9-periods, whose elements read 

Aat = f dut, Bst = f dut, s,t ^ 1, . . . ,{n - l)m , (C.3) 

respectively, can be written in terms of the following mxm matrices {k — 1, . . . , 71— 1) 

iAk)ij = f , (6fe)y = * o?Mj+™(fc_i) , ij = l,...,m, (C.4) 

as follows 

A = diag(y^i, . . .,An~i)'R-A , B = diag(Si, . . . ,6„_l)7^e , 

where diag(^i, . . . ,An-i) and diag(Ki, . . . ,Bn-i) are block diagonal. The matrices 
TZj, and TZjs are defined as TZ^ = TZji id^ and TZjg = TZts (E) id^ (id™ is the mxm 
identity matrix), where the (n — 1) x (n — 1) matrices TZj[ and TZg have the following 
elements 

(^^)fc, EE /(I-'-) , (^g),,^^i^(p-fc'-_i), p^e^. (C.5) 

1 — p*^ 

Introducing the normalised holomorphic differentials dv — {dvi, . . . , dv(^n-i)m) = 
duA~^, allows to define the Riemann period matrix as 



11^4= * dvt, s,t = 1, . . . , (n - 1)™, (C.6) 

which turns out to be 

n = diag(^r''8i, . . .,A-\Bn-i)nB . (C.7) 

Given these definitions, in [36 the following Thomae-type formula (see (5.28) of [36] ) 
has been proven 

n (---..)^^""" n (^2...-.2..o^^"""-(c.8) 



l<.i<j<m 0<i<j<.m 
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In our case, we have m — I i.e. four points 21,22,23, and 24; thus p{z) — 
(2 — 2i)(2 — 23) and q{z) — {z ~ Z2){z — 24). In the rhs of (|C.8|) . the product of 
{z2i — Z2j) is 1 because 1 < i < j < 1 cannot be fulfiUed, while the product of 
{z2i+i — Z2j+i) has only one term, which is 1. The two independent cycles ai and /3i 
coincide respectively with Ci and C2 in Fig. lAll The matrices Ak are one-by-one, and 
therefore 

n = -R-j^ diag(6iMi, . . . , Bn~i/An-i) TZb ■ 

Moreover, TZj\ — TZj[ and TZ^ = TZ^ and (|C.2[) reduces to duk{z) — ojk{z)dz. Thus, we 
have that Ak and Bk are given by the integrals in Eqs. (|A.22p . We have the following 
expression for 11 

n = Tt;^! diag(. . . , p'^/^pk/n , . . .) ■ (C.9) 
By introducing the (71 — 1) x (71 — 1) matrix Mi having all the elements equal to 1, we 
can write 7t^^ and TZb in terms of the matrix U as follows 

TIa = ^ + Ml) U diag(. . . , p-^ . . .) , (C.IO) 

Ub =\/Hdiag(^...,^^,...^ C/t(id„_i+Mi) . (C.ll) 

Now we observe that the matrix (|C.9p is related to T of Eq. ^ as 

n = A-^TA~\ (C.12) 

where A = 2U'^HU with H = diag(. . . , sin(7rfc/n), . . .). As for the rhs of ([Oil , it 
reduces to (notice that ]\l'Zl p^'^ = {^i) ) 

n ^[\{2F,{k/nA-klnA;x)\ . (C.13) 

k = l ^ ^ \fc = l / 

Using finally Q{0\A-^ TA-^) = e(0|r) Eq. ([RS]) . we have therefore proved ([6T|). 

Appendix D. Three analytic continuations more 

Appendix D.l. The analytic case x = 1/2 

For X = 1/2, the hypergeometric function can be written in terms of T function as 

Fk/n{^/2) = _ fc/(2n))r(fc/(2n) + 1/2) ' 
so that the Z?„ simplifies to 

Ai(l/2) = logTT ^'fe/n(l/2) = log ; . (D.2) 

H ^ ^nr=i[r(l/2 + V(2n))]2 

We can then employ the following integral representation of the logarithm of the F 
function 

dte^ / e-("-i)* - 1 
~ V 1 -e- 







logF(2)=/ -^( , +^-l), (D.3) 



and so 

'1 I \ r°° die-' /g-(;/(2«)-i/2)t_i , 
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The resulting series is easily summed 
and the derivative wrt n taken 

r.n.l.^ 1, /•°°die-*e3*/2(2t-5) + 5e* + e*/2-l , , 

^^(l/2)--2l^^ + j/ (e^ - ^ 1) • 

The integral is convergent, but the various pieces in which it can be divided are not 
and so a lot of care must be used to perform it. The easiest way we find out to make 
the integral is to make a sort of dimensional regularization of the term — s- t~^^'^ , 
making then the various integrals that are finite for large enough v and then expanding 
the result close to = where all the divergences cancel leaving (after long algebra) 
the finite result 

D[{l/2) = - 7_E + Iutt) ~ -0.216243 , (D.7) 

where is the Euler 7 constant. 

Appendix D.2. Perturbative expansion in x 

The easiest way to perform the analytic continuation is to expand in power of x 
the function i^^/„(a;). From the well-known series of the hypergeometric function the 
coefficient of this expansion at order p is a polynomial of order 2p in k/n. The product 
over k from 1 to n — 1 can be done for very large orders and the derivative for n = 1 
taken. With a little help from Mathematica we obtain the following perturbative series 
for D[(x)B 

x 2 p 26 o 2 4 52 . 302 76 7 398 . 

- D (a;) = - H x^ H x'^ H x^ H x^ H x^ ^ x^ H x^ 

' 3 15 315 35 1155 9009 2145 255255 

327128 „ 18047684 31378136 ^ ^, 
+ 2297295^ - 24249225^ + ^31173^^ + ^ ' ^^'^^ 
It is easy to obtain as many terms of this series as we want, but this is useless if 
we perform a direct summation. In fact, it is evident that the terms up are all 
positive and quickly decreasing in magnitude, giving the impression of a convergent 
series. Unfortunately from the term they start increasing abruptly and oscillating 
signaling that D'^{x) is an asymptotic series. It would be possible to resum this series 
a'la Borel, but this is of no use in view of the exact result we obtained in the main text. 
This result, that is very easy to obtain, has been a fundamental cross-check for other 
relations we found. It also provides a very good estimation of D'i{x) for x < 1/2. For 
example, summing the first eight terms for x — 1/2 we obtain D[{l/2) ^ —0.216102, 
deviating only of the 0.06% from the exact result above. 

I A numerology digression: The denominators in this Taylor series are integer multiple of the ones 
of tan X up to order 10 when D'^ (x) starts oscillating. On the Sloane on-line encyclopedia of integer 
sequences, another "look-a-like the denominators in Taylor series for tana;" can be found, (sequence 
number A156769) that starts differentiating from tanx at the eleventh term (excluding x), exactly 
like this one. 
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Appendix D.3. Analytic continuation for small x 

In this appendix we report the details of the analytic continuation in the small x 
regime. We need to continue to general complex n the sum P„ in Eq. ([55]) . It is 
instructive to consider first two special cases that are easily worked out. For a = \ 
we have 

n-l . , 

/=i [sm(7r//n)] b 

so that P[ = 1/3 (but we remember that this number has no physical meaning, because 
it is exactly canceled by the denominator to give T[{x) = 0). 
The other easily solvable case is a = 1/2, for which we have 

■n-l J I 

sm (7r( / n) 

We can use 

(1 + x)'^ smTT/i 
with /i = l/n, to have 

1 dx 1 /■°° dx x-x^l'^ 



Trio + 71- 7o (1 + x)2 xV" - 1 



(cotha;tanh(nx) — l)c?x , (D.12) 







n 

that is the desired analytic continuation. From this we have 

P'i = - 7 TT -^dx^- — — dx = -.(D.13) 

""Jo (a; — 1)(1 + x)^ TT Jo smhxcoshx 8 

It is also possible to give other simple formulas for all integer values of 2a, but they 
are not values of physical interest. Let us just mention that for a = we trivially 
have P„ = (ti - l)/2 with P{ = 1/2. 

For general < a < 1, the only strategy we found is to expand the argument of 
the sum using 

oo 

T^-^'-'Y.P>^^"'^ (D.14) 
smxP ^-^ 

^ ' k=0 

with pk known (also to Mathematica) . Using this expansion we have 

-,n— 1 ,/ _,oo iin^7 \ n — 1 



- oo 

-E^4- 



l+2(/c — ct) , , , 



where the harmonic number i/r'''' is the natural analytic continuation of the above 
sum (defined e.g. in terms of the Riemann ^ as iji'*'' = C(/i) — C(/i, z + 1)). This sum is 
still not analytically possible, and does not provide a proper analytic continuation yet. 
In fact, it is easy to check that only for integer n it gives a convergent sum (that can 




be simply truncated at a given order to have the desired precision). For non-integral 
values, the sum is asymptotic and must be resummed. We explicitly show how to do 
this only for P{ needed for the entanglement entropy. Taking the derivative wrt to n, 
we formally get 

oc 

P[ = ^pfc7r2('=-")(l - 2(fc - «))C(2(a - k)) . (D.16) 

fc=0 

This sum is asymptotic. We then introduce the function 

oo oc 

V{z) = ^pfcz2(fe-a)(^ _ 2(fc - a))C(2(a - k)) ^ ^^^^^(/c-a) ^ p^^^) 

fc=0 fc=0 

with the coefficient Pk growing like a factorial for large k. By definition we have 
P[ = 'P{n). The a-Borel transform is 

oo p 

that provides a convergent sum. In case we would have been able to perform this sum 
analytically, the original function would be given by the anti-Borel transform 

V{z) = - e-'^H-^"B-p{t)dt, (D.19) 
z Jo 

as it can be checked by expanding in z. (We used the a-Borel transform to cancel 
the effect of the singularity of the Borel transform in the origin, however this has 
no importance for the numerical results). The analytic sum can not be performed, 
and so one should find an effective approximation of the sum (jP.lSp that makes the 
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integral (|D.19p finite. There are several standard methods to do this and we use the 
Pade approximation of the series, that is a ratio of polynomial of order N and D 
for numerator and denominator respectively. We consider in the integral (|D.19p large 
enough N and D, so that the value of P{ does not change (at the required accuracy) 
when still increasing them. Within this procedure we got the values of P{ for any 
< a < 1 that are reported in the Fig. IDll together with Pn/{n — 1) for integer n. 
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